Metastable Intermediates in the Condensation of Semiflexible Polymers 
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Motivated by results from an earlier Brownian Dynamics (BD) simulation for the collapse of a 
single, stiff polymer in a poor solvent [B. Schnurr, F. C. MacKintosh, and D. R. M. Williams, 
Europhys. Lett. 51(3), 279 (2000)] we calculate the conformational energies of the intermediate 
(racquet) states suggested by the simulations. In the absence of thermal fluctuations (at zero 
temperature) the annealed shapes of these intermediates are well-defined in certain limits, with 
their major structural elements given by a particular case of Euler's elastica. In appropriate units, 
a diagram emerges which displays the relative stability of all states, tori and racquets. We conclude 
that, in marked contrast to the collapse of flexible polymers, the condensation of semiflexible or 
stiff polymers generically proceeds via a cascade through metastable intermediates, the racquets, 
towards a ground state, the torus or ring, as seen in the dynamical simulations. 

PACS numbers: 87.15.He, 36.20.Ey, 87.15.-v 



I. INTRODUCTION 

The conformation of individual polymer chains de- 
pends on the properties of their environment, i.e. the 
solvent In the presence of a poor solvent, isolated 

polymer chains tend to collapse toward compact states, 
in which polymer-solvent contacts are minimized. For 
flexible polymers, the kinetics of this coil-globule tran- 
sition have been the subject of much research over the 
past few decades (|, ||, [| @, ||, [J. The kinetic path- 
way for flexible polymer collapse has only recently been 
experimentally confirmed to involve the formation of a 
pearl necklace and the gradual diffusion of large pearls 
from the chain ends |L0| pd| . 

In contrast to the flexible case, many polymers exhibit 
substantial bending stiffness, thus adding the (oppos- 
ing) tendency to form extended structures. This makes 
a compact globule energetically unfavorable for semi- 
flexible polymers because compact globules involve large 
amounts of bending. Such chains are described by the 
persistent or worm-like chain (WLC) model pj, exam- 
ples of which include predominantly biopolymers {e.g. 
F-actin and DNA) but also some synthetic polymers 
(e.g. kevlar). The balance between the tendencies to 
straighten the filament (due to a bending energy) and 
to condense it (due to an effective short-range attraction 
or poor solvent) is at the heart of the condensation of 
semiflexible polymers. 

The apparent equilibrium collapsed state for semiflexi- 
ble polymers is well-known: chains with significant bend- 



ing stiffness can form rings or toroids to avoid incurring 
the large bending penalty of a spherical shell or a glob- 
ule. This condensed state has been suggested and studied 



theoretically 
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demonstrated in a variety 



of experimental systems ]l6| [r7[ [Lsi [l9|, |2fX_ 2l| , and con- 
firmed by computer simulation [p2[ |23|, ]24, g5[. Theo- 
retical work has predominantly addressed structural fea- 
tures such as the detailed packing of filaments J2G|, |27|, 
p8| , |29| , |30| , , while dynamical simulations and Atomic 
Force Microscopy (AFM) among others have increasingly 
focused on kinetic aspects and condensation intermedi- 
ates 

A particular set of recent dynamical simulations of iso- 
lated chains Q has strongly suggested a possible (and 
in fact generic) pathway for the collapse of semiflexible 
polymers. These simulations showed not only the even- 
tual formation of tori from extended chains quenched in 
poor solvent but demonstrated a series of long-lived, par- 
tially collapsed intermediate states. Very similar chain 
morphologies (our racquet states) also appear in other 
simulation work |24| and AFM studies of DNA condensa- 
tion JT^j. Motivated by these results, we develop and ana- 
lyze a hierarchical family of metastable racquet states. In 
particular, we demonstrate that their relative conforma- 
tional energies are consistent with the role they play in 
the simulations: they form an energetically driven cas- 
cade of increasingly compact conformations with sharp 
transitions between them. 

We begin in Section II with a brief summary of the 
dynamical simulation results |34|| which motivated this 
analysis. Section III addresses the morphology and evo- 
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lution of the shapes to be analyzed in the remainder of 
the paper. Our approach to calculating the surface con- 
tributions to the conformational energies is developed in 
Section [V , followed by the two main sections containing 



a detailed analysis of torus and racquet states (Sections 



2 



and VI respectively). Section VII finally compares their 
relative stability and discusses the qualitative agreement 
with the dynamical simulation results we set out to un- 
derstand. 



II. BRIEF REVIEW OF THE DYNAMICAL 
SIMULATION RESULTS 

The work described in Ref. [|4| applied a standard 
Brownian Dynamics (BD) algorithm B to a bead-and- 
spring model of a single polymer chain in the plane to 
capture the most general features of a rather complex 
and biologically important process, the condensation of 
DNA. The technical details of that study are discussed 
elsewhere |34|, [35j| . Here, we merely sketch the gross fea- 
tures and the generic results that motivated our work in 
this paper. 

The dynamical evolution of a simulated chain followed 
a Langevin equation of the form 



_dxi dll , . ^ 



(1) 



for each bead i, where £ is the coefficient of viscous drag 
(Fvisc = £v) and r\ the random noise. Each bead is 
displaced by Axi — (F x J^)At during a time step At. 
The potential U contains all interactions internal to the 
chain, including the bending energy, a short-range attrac- 
tive interaction between beads (mimicking poor solvent 
conditions) and a very stiff longitudinal compliance. Af- 
ter thermalization of each chain, the solvent quality was 
quenched at t = 0. 

Previous work p4| , |35| showed the typical dynam- 
ical evolution of a relatively short chain (a few per- 
sistence lengths) as a progression through well-defined 
stages identified by three types of conformations: ex- 
tended chain with thermal undulations, various racquet 
states (see Figs. |l| and ^), and the torus or ring. We also 
pointed out that the end-to-end distance of the filament 
as a function of time changes sharply with the confor- 
mational transitions between states. It is important to 
note that the described conformations persist in time, as 
seen by quasi-plateaus in the end-to-end distance evo- 
lution, each lasting for the considerable time of about 
10 6 BD steps, about one tenth of the entire condensation 
event. We can roughly estimate the correspondence be- 
tween simulation steps and physical time for a particular 
system. To do this, we express the link length as a frac- 
tion of the persistence length and substitute for the local 
drag coefficient, assuming the viscosity of water. For F- 
actin, such an estimate suggests that an entire simulation 
of 10 7 BD steps models a filament for about 0.1 seconds. 
For DNA, this interval corresponds to a fraction of a mil- 
lisecond. 

Temporal persistence of racquet structures was seen 
throughout the simulations, suggesting that metastable 
intermediates are a general feature in this collapse. Pre- 
sumably, energy barriers between intermediates are re- 



sponsible for their local (meta-)stability but we have not 
attempted to estimate their size. 



III. IDENTIFICATION OF INTERMEDIATE 
STATES 

The dynamical simulations [Q suggest that the mech- 
anism of collapse of semiflexible chains generically in- 
volves transitions through a series of long-lived interme- 
diate states. In the absence of thermal fluctuations these 
intermediates anneal to certain underlying shapes which 
are well-defined and allow a straight-forward calculation 
of their conformational energies. The crucial element in 
the underlying shapes is a characteristic looped section 
which we call a racquet head. For the single racquet, the 
shape of the head (see Fig. [I]) was produced in the simu- 
lation by annealing. While missing the effects of thermal 
undulations, our calculations of the conformational en- 
ergies and detailed shapes of the annealed intermediates 
provide an insightful framework for understanding the 
simulation results. 




FIG. 1: Annealed shape of a racquet head from the BD simu- 
lation, achieved by slowly lowering the effective temperature 
once the structure has formed. The shape coincides with the 
analytical curve to within a line width. 

In order to simplify the reference to specific states, 
we label all racquet states by their number of looped 
sections. Thus the rod is the N = state, the racquet 
with a single loop at one end the N = 1, and so forth, as 
indicated in Fig. |[ We refer to the loop formed at the 
ends of the structure as the "head" and to the bundle 
of filaments connecting heads as the "neck" . For the 
moment we neglect the more subtle question of the exact 
location of filament ends. Naively one might assume that 
filament ends coincide with the ends of the neck, since it 
is straight; that is, the ends of the filament are expected 
to span the entire neck as they incur no bending penalty, 
but generally gain from increasing their overlap. 

The picture as described provides an adequate starting 
point for labeling the states we consider here. Among the 
shapes indicated in Fig. || we distinguish two basic rac- 
quet symmetries: even and odd total numbers of heads 
N. For racquets with even N, the number of overlaps in 
the head sections is equal on the two sides, p — q = N/2. 
For racquets with odd TV, one side (we arbitrarily call it 
the left, following Fig. 0) has one less p — (N — l)/2 than 
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the other: q = (N+l)/2. As a consequence, the filament 
ends of an even racquet are on opposite sides of the neck. 

Note that the dynamical simulations modeled the case 
of fixed experimental conditions after the solvent quench. 
The polymer chain is merely exploring a given conforma- 
tional energy landscape via thermal fluctuations. How- 
ever, it can be instructive to consider as a Gedanken- 
experiment the case of variable filament length, and we 
use this perspective in our discussion. A (reduced) chain 
length is also a natural independent variable for the pre- 
sentation and comparison of states. In this alternative 
perspective, the evolution of shapes starts with a short 
filament that gradually lengthens. At first, only the neck 
grows until the formation of a new head is favored. The 
incremental unit of growth between conformations thus 
consists of one head plus one neck segment. This pro- 
cedure can be conti nued to arbitrary N given enough 
filament. In Section IV C| we will see that the appropri- 
ate formulation of the problem accomplishes changes in 
the (effective) filament length by adjusting the solvent 
quality instead of the actual chain length. 



simulations |34j confirm this. Most of the chains stud- 
ied there were relatively short (less than 10 persistence 
lengths) though a few examples (see Fig. ||) of longer 
chains showed qualitatively similar behavior but with in- 
creased complexity, such as the display of superstructures 
of racquets within racquets. 




FIG. 3: Early stages in the evolution of a long chain (300-mer, 
roughly 20 £ p ) showing combinations of the conformational 
elements seen in shorter chains. 
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FIG. 2: Schematic "family" of racquet states. The rod can 
be thought of as a trivial racquet without head (N = 0). 
All subsequent states are labeled by their numbers of head 
sections. 

In the simulations, the actual transition into the torus 
state could not be resolved in detail. It is clear however 
that there are in principle at least two ways in which a 
loop can form: two chain segments can meet with their 
tangent angles at an angle of 7r or 2tt. The former case 
leads to a racquet head, while the latter makes a ring 
that allows the chain to wind up into a torus directly. 
Since it is more likely for a stiff chain to bend into the 
smaller angle, one would expect the transition to racquets 
to be favored, at least for short chains. Statistically, the 



IV. CONFORMATIONAL ENERGIES 

Having discussed the basic morphology of the inter- 
mediate states, we turn to the calculation of their con- 
formational energies in the absence of thermal undula- 
tions. Racquet and torus conformations at zero temper- 
ature can be thought of as the underlying shapes, which 
become modified by fluctuations at finite temperatures. 
Apart from the bending contributions to the conforma- 
tional energies, we need to describe the nature of the 
surface energy by which we model the poor solvent con- 
ditions that induce condensation. This surface energy 
assumes the packing or bundling of filaments in a hexag- 
onal lattice (in cross section) and distinguishes between 
polymer and solvent exposure. The local arrangement 
into a hexagonal columnar phase has been confirmed for 
example by x-ray diffraction applied to bundles of DNA 
and other charged polymer chains |37], [58], EJ 
and the detailed structure within bundles of semiflexible 
polymer chains has been studied theoretically |4l], Q . 

Our calculations describe the simplified model with fil- 
aments of vanishing thickness (compared to other length 
scales in the problem), while we model their packing on 
a perfect hexagonal lattice. Thus we do not take into ac- 
count any winding defects due to topological constraints 
or variations in curvature due to the finite filament thick- 
ness. As implicit in the description of a worm-like chain, 
we assume a uniform bending modulus. In the torus 
state, such an ideal chain forms a circular ring with a sin- 
gle radius of curvature. Furthermore, we assume for both 
racquet heads and torus that partial filament overhang, 
an effective non-uniformity in the bending modulus, does 
not change their shape but only their size. Neglecting 
these higher order corrections is certainly justified in the 
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limit that the bundling number N gets large. 



A. Surface Energy of Hexagonal Bundles 

In a hexagonally close-packed bundle, each filament in 
cross section can be thought of as having six sites occu- 
pied by either solvent or polymer. The poor solvent low- 
ers the energy for polymer-polymer relative to polymer- 
solvent contacts. To express the fact that there is a rel- 
ative energetic advantage for filaments to bundle versus 
being exposed to solvent, we explicitly evaluate the total 
number of solvent-exposed sites and express the energy 
as a surface tension. 

Particular surface energies are evaluated as follows. To 
find the coordination number oln for an N bundle, con- 
sider the total number of surfaces or binding sites in the 
bundle with hexagonal order (6N). This number is pro- 
portional to the energy of N individual filaments com- 
pletely exposed to solvent. To account for the effect of 
bundling, we note that a bond corresponds to the merg- 
ing of two binding sites on neighboring filaments. We 
thus subtract the number of bonds formed from only half 
the number of sites (3JV) to find the coordination num- 
ber oln- As an example of this numerology consider the 
cases for N = 5 and 10: for 5 filaments there are 7 bonds, 
resulting in a coordination of 8, while 10 filaments make 
19 bonds and thus have a coordination of 11. Multiply- 
ing by the surface tension parameter 7 finally yields the 
surface energy per unit length for such a bundle. 



B. Filled Shells and "Magic" Numbers 

Differences between subsequent ajy are always cither 
zero or one (except between N = 1 and 2). This creates 
non-uniformities in the effective binding strengths per 
unit length, thus favoring particular bundling numbers. 
We expect this effect for N with the same coordination 
as their predecessor (ajv = ajv-i) but a coordination 
of one less than the (N + 1) bundle («n + 1 = aw+i). 
This is the case whenever an added filament adds three 
instead of two bonds, thereby filling a shell. Examples of 
this situation are found for N = 7, 10, 12, 14, 16, 19, 21 . . . 
and we refer to them as "magic" numbers or filled shells. 
Cross sections of magic number bundles correspond to 
arrangements with high degrees of symmetry, as shown 
in Fig. [|. 
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A few "supermagic" numbers (N — 7, 19, 37 ... ) rep- 
resent bundles with the special symmetry of the equilat- 
eral hexagon; we will not treat these cases separately. In 
agreement with Pereira and Williams (43) we find that 
bundles of magic numbers (and particularly the super- 
magic ones) play the role of preferred states with in- 
creased stability. 



C. "Condensation" Length and Energy 

The formulation of the problem as presented con- 
tains a characteristic length scale which greatly simplifies 
the discussion and presentation of our results. Balanc- 
ing expressions for typical bending and surface energies 
(k/L jL) for a given filament length L defines a mea- 
sure we call the condensation length L c = \J fc/7. Its role 
in the behavior of a chain under particular conditions is 
the following: given the physical parameters k and 7 a fil- 
ament much shorter than L c will rarely self-intersect and 
therefore typically form an extended structure, while one 
much longer than L c is likely dominated by overlaps and 
will form collapsed or at least partially collapsed (inter- 
mediate) structures. 

Another combination of the two basic parameters k 
and 7 sets an analogous energy scale, the condensation 
energy U c = ^/kt/- With these measures, all confor- 
mational energies Un can be presented in dimensionless 
units, where physical energies and lengths are normal- 
ized by their condensation values: u/v = Un/U c and 
A = L/L c . This formulation also provides a conve- 
nient (experimental) realization of "changing the filament 
length" . We can vary the reduced length A by adjusting 
the values of k and 7 independently. 



V. TORUS STATES 

We describe a torus by the following two (dimension- 
less) variables: a filament of length A is wound into a cir- 
cle of constant radius p, as shown in Fig. |j| In general, 
the torus can have any number N of complete windings 
(through an angle 2tt) and an amount of extra overhang 
a subject to the condition a < 2irp. Since the entire 
filament contour length A has a constant curvature, the 
bending contribution to the conformational energy is al- 
ways A/2,0 2 . 




FIG. 4: Bundle cross sections of the lowest magic num- 
bers on a hexagonal lattice. Note the arrangements of 
the "supermagic" numbers 7 and 19 into perfect hexagons. 



FIG. 5: Sketch of a generic torus (here a 1+ in our labeling 
scheme) of radius p and overhang a. 
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We find it convenient to distinguish forms of the torus 
with different numbers of complete revolution N defined 
as the largest integer in X/2np. Any non- integer por- 
tion of this ratio represents overhang of filament beyond 
complete windings, defined as a = A — N(2np). Our dis- 
tinction of different tori by N naturally separates cases 
with different coordination numbers and thus different 
surface energies. In anticipation of an important distinc- 
tion that emerges, we call a torus without extra overhang 
(<t = 0) an "exact N" while we refer to the generic torus 
with finite overhang (a ^ 0) as an U N+" . 

For the torus as described, we can then write down the 
following expression for its total conformational energy 

i$ rus (A, P, o") = ^2 + 277 P + a ^ aN + 1 ~ aN ^ ( 2 ) 

with the common bending term followed by two surface 
terms describing the contributions from the complete N- 
fold ring and the extra piece of overhang tr, respectively. 

Substituting for a leads to the torus energy in terms of 
A and p only, which allows us to find the equilibrium size 
or radius Pn{X) for a particular state N by minimiza- 
tion with respect to p: note that d 2 u t § TUS /dp 2 = 3X/p 4 
is positive everywhere. Resubstitution of pn(X) yields 
the conformational energy u^ rus (A) in terms of the sin- 
gle variable A. The expressions found in this way are 
valid in the ranges of A between N and N + 1 times 
the circumference lixp^. However, a real solution for 
this equilibrium size need not exist. This happens ex- 
actly for the magic numbers with N > 12. When a 
real solution for the equilibrium radius does exist, the re- 
sulting energy has two terms: one proportional to A 1 / 3 , 
the other to A. The coefficient of the linear term is 
the combination of coordination numbers (a^+i — ajv) 
which can vanish for N just below the magic numbers 
(at N — 6,9,11,13,...) leaving these cases with the 
functional dependence A 1 / 3 only. Another consequence 
of the numerology of hexagonal packing is that differ- 
ent states N can share the same energy expressions. 
Examples are the series N — (2,3,4,5) and the pairs 
N = (7, 8), (17, 18), (22, 23), (25, 26), (28, 29)... These 
cases form a particular class of transitions where a grows 
continuously with A. 

A. Stability 

The above method for finding the optimal torus sizes 
by minimizing the conformational energy represents the 
conventional approach to determine metastability. How- 
ever, there is a somewhat unusual aspect to the problem 
at hand. The energy expressions for torus states with 
different N are in general not the same. This introduces 
discontinuities in the form of the energy between adjacent 
states. Consequently, there are not only the conventional 
minima identified by their vanishing slopes but also an- 
other class of solutions with discontinuities in slope at 
points where the energy expressions to the left and right 



differ due to the filament coordination. These are not 
minima in the usual sense (for instance, they are not 
locally quadratic minima); they are stabilized by finite 
slopes on both sides and do not have the usual signature 
of a vanishing slope. Our results for the tori are displayed 
in Fig. |. 




20 40 60 80 100 
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FIG. 6: Conformational energies of the torus states as a func- 
tion of filament length in reduced units. The thinner lines 
indicate (metastable) solutions in regions where they are not 
the ground state, which in turn is indicated by bold segments. 
The dashed line shows the large N solution calculated in Sec- 
tion |vc| 

In order to establish the metastability of the tori of dif- 
ferent N in more detail, we consider the behavior of the 
energy derivatives with respect to the radius du*§ lus /dp 
evaluated at the radii where the exact N and N + 1 form. 
These derivatives are all monotonic functions (functional 
dependence: —A -2 ) with at most a single zero indicating 
a limit of metastability. Around these zeroes, the deriva- 
tive is generically negative to the left and positive to the 
right. When negative, the energy is lowered by increasing 
the radius p, thus driving any overhang a to vanish and 
making the torus an exact N. When positive, the oppo- 
site is true, driving a to grow, making the torus an N+. 
Note that there are cases (notably again for the magic 
numbers with N > 12) where duN/dp is negative ev- 
erywhere. These cases form an important class in which 
tori never evolve (with increasing A) towards states with 
finite overhangs: they remain metastable (with complete 
or exact overlap) for all lengths beyond some lower limit. 
The relative positions of the zeroes in the energy deriva- 
tives combine in two fundamental ways, resulting in ex- 
act and N+ tori for various ranges of A. A more detailed 
discussion of the various cases can be found in |35[ . 

The rod (N — 0) is of course a special (trivial) case 
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without any bending contribution. Due to the absence 
of any competition between bending and self-affinity, the 
rod is (at least) metastable for all lengths and follows a 
straight line. The lower limit of the 1+ state is given 
by the circumference of a single ring that just closes 
(A = 27r), a circle with the radius of a condensation 
length. The subsequent small N states show variations 
depending on the numerology of the hexagonal coordi- 
nation numbers. For larger N a perhaps generic type of 
series emerges where N+ become exact N + 1 which re- 
main metastable to infinity. Thus, in contrast to what 
we have emphasized here by treating only cases with rel- 
atively small N in detail, there appear to be only the few 
tori with N below 12 that show variations to the generic 
pattern of exact (magic N) tori without extra overhangs. 

A direct comparison between relevant branches of the 
solutions over regions of A provides the transition points 
between stable (or ground) states, as summarized in Ta- 
ble | 



TABLE I: Transition points for the lowest energy states 
(ground states) up to N = 24. Only for the "shortest" chains 
(A up to 11.543) is the rod stable to the torus. 



State Labels 


Transition Points 


0+ . 1 + 


11.543 


1+ > 2+ 


12.957 


2+ . 3+ 


18.850 


3+ > 4+ 


29.021 


4H > 7 


38.871 


7 > 10 


73.625 


10 > 12 


93.195 


12 > 14 


119.876 


14 » 16 


148.687 


16 > 19 


155.672 


19 > 24 


228.700 



It appears that the majority of stable torus states (per- 
haps all TV > 12) are exact states over their entire range. 
Their labels N are a subset of the magic numbers. 



B. Discreteness 

Another notable result are the discontinuities due to 
the hexagonal packing and its discrete coordination num- 
bers. We might have expected the overhang a and the 
torus size p to be continuous with changes in A. At least 
for small N we find instead that small changes in A can 
cause discrete jumps in the size of the ring p. This char- 
acteristic has previously been described by Pereira and 
Williams (43) . To what extent these effects are experi- 
mentally observable is not known. As N grows large, the 
effect should weaken and ultimately disappear altogether. 

Discrete jumps are perhaps most prominently dis- 
played as discontinuities in the torus size. Fig. ^ shows 
the sizes or radii pjv(A) of the ground states as a function 
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FIG. 7: Radii of the stable torus states (bold) as a function of 
filament length in reduced units, showing discrete transitions 
between the various series of states. The series shown are 
labeled by the states at their lower extremes. For comparison, 
the continuous large N solution is shown as a dashed line. 



of A. The first bold curve segment starts where the fila- 
ment first makes a stable 1+. Note that the functional 
dependence of the first two segments is different from the 
subsequent (linear) ones. The first two series evolve con- 
tinuously according to their equilibrium solution for p^ 
with a functional dependence of A 1 / 3 . Their prefactors 
are determined by some combination of the appropriate 
coordination numbers. By contrast, all subsequent seg- 
ments are due to solutions which are constrained to be 
exact (by virtue of the magic numbers) and thus have the 
linear dependence of \/2irN . The length of the various 
segments indicates the stability of the states they repre- 
sent. Clearly, states with supermagic bundling numbers 
(the figure only shows N — 7 and 19) are especially sta- 
ble. The dashed line indicates the (continuous) solution 
found in the limit of large N, as discussed in the following 
section. 



C. Large N Limit and Scaling 

The scaling argument of the torus size with filament 
length goes back to Ubbink and Odijk j pg} ]. We sketch 
a similar argument here in order to compare it in Sec- 
tion VI D with the analogous argument for the racquets. 
First we give the straightforward scaling argument. In a 
second pass, we then determine the prefactors based on 
the more accurate, hexagonally faceted cross section of 
the torus. 
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In the limit as TV grows large, we can neglect such de- 
tails as partial overhangs (finite a) since differences be- 
tween TV and TV + 1 vanish as 1/TV. We assume first that 
the torus is a perfect cylinder with circular cross section. 
It grows as TV, the number of filaments wound around its 
circumference. Thus, we expect the total torus surface 
area to scale as p\J~N . Substituting for the radius in a 
scaling sense (p ~ A/iV) we find that the conformational 
energy has two terms: one proportional to A/TV 1 / 2 , the 
other to N 2 /X. Minimization with respect to TV (im- 
plicitly letting TV be a continuous variable) yields the 
following set of scaling relations. 



prefactors are then 



TV 



Po 



A 4/5 
A l/5 
A 3/5 



(3a) 
(3b) 
(3c) 



To find the prefactors, we need to consider a geometri- 
cally more careful treatment. We assume that the torus 
formed has both a perfect hexagonal cross section and an 
integer winding number (no partial overhangs). These 
assumptions are reasonable: it can be shown by direct 
calculation that the surface tension for a fixed number of 
filaments on a triangular lattice is smaller for the hexag- 
onal than for the circular cross section. This is analogous 
to a Wulff construction [Q which captures, for instance, 
the faceting of crystals in solid state physics. The inte- 
ger winding number is justified since the difference in the 
surface energies between TV and TV + 1 filaments vanishes 
as TV becomes very large. 

For an (equilaterally) hexagonal cross section with its 
symmetries, we can determine the following relationships 
geometrically. As a characteristic for the size of the 
hexagon, we label the integer number of lattice spacings 
on a side by m. The counting of lattice sites (or filaments) 
in such a hexagonal bundle is TV = 3m 2 to leading order 
in m. Proper counting adds a linear and a constant term: 
TV = 3m 2 + 3m + 1 but in the limit of large TV we keep 
only the leading order in m. 

We find the surface energy of such a bundle by counting 
solvent-exposed filament sites. A filament at an edge 
(of which there are m — 1) exposes 2 sites, while one 
at a corner exposes 3. Taking into account the 6-fold 
symmetry of the hexagon, this results in 12m + 6 exposed 
sites on the surface. Substitution then yields the limiting 
coordination number (a surface energy per unit length) 
ckoo = 2v / 3~TV for a hexagonal bundle of TV filaments. 
This coordination number also provides the prefactors 
for the conformational energy of the torus in the limit of 
large TV. 

2V3A 2tt 2 TV 2 

= (4) 

The expressions analogous to Eqs. || with geometrical 



TV = 



31/5 



A 4/5 « 0.286 A 4/5 



(5a) 



(2tt) 4 /5 

Poo = (6tt)- 1/5 A 1 / 5 « 0.556 A 1/5 (5b) 
A 3/5 « 8.092 A 3/5 . (5c) 



5(3^) 2 / 5 

2 3 / 5 



The last two expressions are shown in Figs. ^ and |?] as 
dashed lines. We see outstanding agreement between the 
large TV limit and the exact solutions down to the lowest 
TV in Fig. |. 



VI. RACQUET STATES 

The racquet conformational energies are made up of 
bending contributions from each of the heads, and sur- 
face contributions from the heads as well as the neck 
region in between. As shown in Fig. ||, the racquets di- 
vide naturally into two groups: those with even and odd 
numbers of heads. 



N+l 



Even 
(N=2) 



11 = N/2 



N+l 



q = p+i 




Odd 

(N=3) 



FIG. 8: Comparison of the structure and labels of generic even 
and odd racquets, represented here by the TV = 2 and 3. All 
the heads on one side are identical; the schematic separates 
head and neck filaments only to indicate their multiplicity. 

In the even case, the number of heads on each side 
equals n = TV/2 by symmetry. In the odd case, we have 
p = (TV — l)/2 heads on the "left" and q = (TV + l)/2 
on the "right" . The labels "left" and "right" are our 
arbitrary naming convention (see Fig. |8). The variables 
p and q for the bundling numbers of the heads always 
differ by 1 (g = p + 1) and sum to TV. The bundling 
number of filaments forming the neck is always TV + 1. 

Given these bundling numbers, the remaining variables 
(in dimensionless units) for the generic racquet are the 
overall filament length A, and the head sizes on the two 
sides (namely the contour lengths of the heads, labeled 
X P and Xg)- So far we have described the racquets with 
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their filament ends coinciding with the ends of the neck. 
However, in general (and in analogy with the overhang 
a in the torus case) we need to allow for the extension 
of these ends into the heads, or the retraction into the 
neck. Lengths of overhang are labeled a p and a q and 
their sign indicates whether they extend into or retract 
back from the heads. The length of the neck r is not an 
independent variable once all other parameters are fixed, 
since the total filament length imposes a constraint. 

For the even racquet with a given A, the number of 
variables reduces to only two. Since the left and right 
heads for even racquets are identical by symmetry, we 
collapse their labels and are left with only one head 
size (xn = Xp ~ Xq) an d a single overhang variable 
(a n = a p — a q ). The overall filament length for the even 
racquet is distributed into A^ on = Nx n + 2ct„ + (N + l)r 
where the terms are ordered as heads, overhang, and 
neck. For the odd racquet, we leave the left and right 
head sizes separate, but require that any overhang be 
symmetrically distributed on the left side. This is not 
the only possible metastable solution, but the one we 
describe here generically as the most symmetrical; we 
discuss the d etails of other possible solutions further in 
Section VIC. In the odd case, the overall chain length 
divides itself into A^ d = px P + qx q + 2<r p + (N + l)r 
where we use the single overhang variable a p to indicate 
that the two possible pieces of overhang are always on 
the left side (see Fig. |). 



400 




100 150 200 
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300 



FIG. 9: Racquet state energies, shown as a function of fila- 
ment length in dimensionless units, form a dense spectrum of 
solutions, each increasing linearly. For comparison, the scal- 
ing solution with proper prefactor in the limit of large N is 
superimposed (dashed line). Note the nearly perfect agree- 
ment of the scaling solution with the lower envelope of the 
racquet states, down to the very lowest values of N. 



By way of a preview, we state here that the racquet 
solutions (see Fig. ^|) differ fundamentally from those of 
the tori. While the size of the torus was found to increase 
as a function of A (up to discontinuities), the sizes of 
the racquet heads (as well as any lengths of overlap) are 
fixed for each state by the local force balance between the 
bundles of filaments making up head and neck. Having 
determined the head sizes and overhangs for a particular 
state, its lower limit of validity Ai ow is found by adding 
up the head sizes and overhangs in the absence of any 
neck at all (r — * 0). This is the minimal filament length 
required to form a particular racquet. For all lengths 
beyond, the racquets remain metastable as their energies 
increase linearly with slope ajy+i/{N + 1). Adding extra 
filament to any racquet configuration only lengthens its 
neck, while the head sizes and any overhang remain fixed. 
As a consequence, all racquets are (at least) metastable 
solutions for any A beyond their lower cutoff Ai ow - What 
remains to formulate the total conformational energies of 
the racquets is the bending contribution due to partial 
overhangs. 



A. Head Shape — an Elastica 

Having identified the racquet head as the distinguish- 
ing common element among the intermediate states, we 
calculate its geometrical shape (see Fig. |l0|) in the ab- 
sence of thermal fluctuations from the bending of a slen- 
der, elastic rod. The expression for this head shape is 
necessary for the determination of bending energies for 
the racquets. The general class of shapes resulting from 
the bending of a slender rod by forces and couples ap- 
plied at its ends only are known as elastica. Such solu- 
tions were first studied by Euler in 1744. The particular 
solution we seek is schematically drawn in Fig. 51 of the 
treatise by Love fl5f . 

To solve for the shape of a racquet head of total contour 
length H we consider the geometry as shown and labeled 
in Fig. [n]. In this section we use physical variables in- 
stead of the dimensionless units introduced previously, 
as they are more intuitive here and allow for dimensional 
analysis. Given the obvious symmetry about the y axis, 
it is sufficient to solve for one half of the racquet head 
only. The tangent angle along the curve increases from 
9 = at the origin O (s = 0), via a maximum at the 
inflection point /, to 9 = ir/2 at Vi (s = H/2, where the 
head joins the neck). Note that there are two points V\ 
and V% at which the tangent is vertical (9 = 7r/2), with 
an inflection point / between them. These three points 
define an additional symmetry (about the inflection point 
/) for the contour between points V\ and Vi. 



Our particular elastica is solved |46 by minimizing 
the WLC (worm-like chain) Hamiltonian subject to the 
boundary condition that the two halves of the head join 
in the neck at x = 0. We impose this constraint by means 
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FIG. 10: Schematic figure of a racquet head with axes appro- 
priate for our calculation. Local tangent angles 6 are mea- 
sured from the x axis. Symmetrical regions along the contour 
s are delimited by solid circles. 

of a Lagrange multiplier £. 

u= r ds {i(^) 2+ccose \ (6) 

Applying Euler's equation to this expression leads to the 
differential equation 

— = -/? 2 sin0( S ) (7) 

where we made the substitution /3 2 = £/k and expressed 
the angle 9(s) explicitly as a function of the contour 
length s. Note that the Lagrange multiplier £ has the 
dimensions of a force and expresses the force required to 
join the two filament bundles in the neck. Eq. ^ can be 
integrated to yield an expression for the curvature along 
the head contour s as a function of the tangent angle 9. 

^ = ^l-*W(0/2) (8) 

Equivalently, one can rearrange terms and express the 
contour length s scaled by f3 in the form of an incomplete 
elliptic integral of the first kind F(4>, k) 

e 

dt 

13s = k = = k F(0/2, k) (9) 

Jo yi-/c 2 sin 2 (t) 

where k is the elliptic modulus (yet to be determined) 
and t an integration variable. This is the parametric so- 
lution of an elastica: it gives the contour length s as a 
function of the tangent angle 0. The expression is mul- 
tivalued over its range, but invertible in certain regions. 



Four such regions are defined by the axial symmetry of 
the head (through the origin and the neck) and the pair 
of inflection points in between. We thus cover the entire 
racquet head in a piecewise fashion, while only two of 
these regions are essentially different: the piece from O 
to V\ and that from V\ to / (with its reflection from / to 
V%). The yet unknown elliptic modulus k for our elastica 
is found from a geometrical constraint. By the symme- 
try between the segments around the inflection point, we 
demand that the x value at point V\ is twice that at the 
inflection point /. Solving the resulting equation numer- 
ically gives the value for the modulus as k = 1.1695. The 
inflection point is identified by the vanishing curvature of 
Eq. H which corresponds to a tangent angle of 9i — 2.052. 

The expression for the curvature (in Eq. ||) allows us 
to evaluate the bending energy of such a racquet head. 
Since the bending energy is an integral over the squared 
curvature, we can use Eq. || to evaluate this energy over 
any segment of the racquet head by integration. This re- 
quires the numerical evaluation of an incomplete elliptic 
integral of the second kind E(<fi, k). 

U=^ fdB \/l-fc 2 sin 2 (0/2) = ^ E(6/2, k) (10) 
k J v k 

Adding up the symmetrical pieces of this solution for the 
entire head yields the total bending energy of a com- 
plete racquet head [7 head = A(k/H) with A representing 
the numerical constant 18.3331. Thus, the bending en- 
ergy in a racquet head depends (apart from the value 
for the bending modulus k) only on its contour length. 
Note also that the bending energy of the racquet head is 
very close to (but slightly below) that of a circular ring 
with the same contour length H (U rm s = 2it 2 k/H w 
19.7392 k/H). Using circular rings as racquet heads 
would thus provide a reasonable approximation for the 
calculation of conformational energies, provided we ne- 
glect the penalty due to the sharp bends at the neck. 

The form of the solution in Eq. || reveals that our rac- 
quet head shape or elastica is unique, in the sense that it 
is independent of any parameters in the problem. Both 
the parametric head shape s(9) and the bending energy 
U(9) are scaled by the factor j3, related to the local force 
balance at V2. The overall size of the resulting shape is 
merely scaled up or down, while its aspect ratio remains. 
Any slender, uniform rod, subject to these boundary con- 
ditions, will assume the described racquet head shape. 
We emphasize that the size of the racquet head does not 
depend on the overall filament length A, unlike in the 
case of the torus. 

This dependence of the head size on the local force bal- 
ance also suggests that the following experiment should 
be possible, at least in principle. Evaluating head sizes 
in a sample of partially condensed filaments would mea- 
sure the local interaction strength between filaments, a 
quantity not easily found by other means. This approach 
assumes of course, that the value for the bending modu- 
lus (or equivalently the persistence length) is known from 
an independent measurement. 
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B. Bending Energy in Racquet Heads 



racquet. 



To evaluate the bending contribution to the conforma- 
tional energies we recall the expression for the bending 
energy in a head of size %. Generalized to an N bun- 
dle (which effectively multiplies the bending modulus k) 
the dimensionless bending energy for an N racquet head 
becomes it^ ad (x) = A(N/x) where A is again the same 
numerical constant evaluated previously from elliptic in- 
tegrals. A stability analysis and numerical minimizations 
found that "perfect" racquets (with a = 0) are the solu- 
tion for only a subset of all racquets. 



X 



FIG. 11: Partial overhang a (solid line) into an existing head 
of size x (dashed). The ratio of a/x defines the fractional 
overhang ?. This schematic shows the typical situation for 
racquet solutions with nonzero overhang, with the filament 
end located somewhere between the inflection and the halfway 
points (solid circles). 



In order to account for partial overhang into the heads, 
we need to generalize the notion of the numerical prefac- 
tor A. This "constant" is really a function of the partial 
overhang. Due to the scale invariance of our elastica, 
it is not surprising that A depends only on the relative 
overhang q = u/x- In terms of <j the four regions are 
delimited by the following values: 0, 0.1627, 0.5, 0.8373, 
and back to 1, measured from the neck. The three in- 
termediate values identify the two inflection points and 
the halfway point (the origin in Fig. [l0|); note that these 
values are measured in the opposite sense from the one 
defined in the figure. Reconstruction of the piecewise 
solutions for any amount of partial overhang a yields 
the expression M partlal (tr, <;) = A(<;)/a with the numerical 
prefactor A = A(<; — 1) generalized to the function A(<;). 



C. Even and Odd Racquets 

The surface energy terms for all racquets consist of 
several terms with different coordination numbers in gen- 
eral. The only term the even and odd cases share is the 
coordination in the neck, whose length r is shared by 
(N + 1) filaments. For the even racquet, symmetry sim- 
plifies the expressions somewhat. In each of its heads, 
we find a length (\n — o~ n ) with coordination a n and the 
overhang piece a n with coordination a„+i while the neck 
has the common coordination chn+i which leads to the 
full expression for the conformational energy of the even 



even 

u N 



A 



N 



Xr. 



+ 2A(a n / Xn ) 



+ 2 [a n (xn - <?n) + a n+1 (a n )} 

+ Ct N+ l(T) 



(11) 



The first two terms are the bending contributions for 
complete heads and partial overhang, while the three fol- 
lowing terms are surface contributions for head segments 
and neck, respectively. For the odd racquets, the expres- 
sion becomes 



odd 



A 



P_ 



q_ 

Xq 



+ 2A(a p /x P ) 



+ MXp - 2<T P ) + a q (x q + 2a q ) 
+ a N+ i(r) . 



(12) 



In the even case, there are only two free variables, Xn 
and a n . We find their optimal values by simultaneous, 
numerical minimization. In the odd case, the situation 
is slightly different, since we lack the symmetry between 
heads. However, we can make use of the fact that both 
filament ends (and therefore any potential overhang a) 
are on the left side. This leaves the right head with 
a well-defined structure of q filaments in the head and 
N + 1 = 2q filaments in the neck. Since the head size is 
solely determined by the respective bundling numbers in 
head and neck, we can determine right head sizes inde- 
pendently of any overhang on the left by minimization in 
terms of the various bundling and coordination numbers 
and A only. 



Xq 



2qA 



2a 



Oi2q 



(13) 



We then numerically minimize over the two remaining 
free variables Xp an d o p of the left head. Plotting con- 
stant energy contours as a function of the two free vari- 
ables generally reveals the approximate location of the 
relevant minimum, and their coordinates were used as a 
starting point for the minimization routine. This pro- 
cedure finds two possible outcomes for both even and 
odd cases. In the simpler case, the energy is minimized 
without overhang (a — 0) and we recover the naively as- 
sumed, perfect racquet structure. In the other case, we 
find a local minimum with respect to a and x f° r finite 
values of overhang. In every case we have checked (up 
to N = 30), these solutions put the fractional overhang 
in the second region of the racquet head, between the 
inflection and halfway points, as indicated in Fig. O. 

Since both filament ends are on the same (left) side for 
the odd racquets, there are several possible configurations 
for overhang to be arranged, as shown in Fig. [l2]. The 
two pieces of equal length a p could be arranged symmet- 
rically on opposite sides of the head. Alternatively, the 
two pieces of overhang can be on the same side, but not 
necessarily of equal length. All the cases we examined 
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FIG. 12: Three possible solutions for partial overhang into 
the left head of an odd racquet. Case (a) is the one described 
in the text. The more general case with different amounts 
of overhang on the same side (b) is always minimized by the 
arrangement in case (c) where the two ends coincide. 



are minimized for one unique value of a, corresponding 
to cases (a) and (c) in Fig. [l^ which turn out to be de- 
generate in energy. In retrospect we were thus justified 
to describe the odd racquet with overhang generically as 
the symmetric case (a), while a second (asymmetrical) 
solution, degenerate in overhang and energy, exists. 

All racquet head sizes Xp an d Xq found either by direct 
calculation or by numerical minimization are displayed 
as a function of the racquet state label N in Fig. [l^ to 
show the general trend and their convergence towards the 
large N solution. Head sizes typically increase with TV, 
though not monotonically, and right heads are typically 
larger than left for the odd racquets. Even racquet heads 
are of the same size, by construction. 
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FIG. 13: Left and right head sizes versus the state label N. 
Pairs of even heads are of the same size, by symmetry. The 
general trend is for heads to grow with TV, if not monotoni- 
cally. Notice the convergence towards the asymptotic solution 
(dashed line, see Eqs. [L4j) with increasing N. 

Since our minimization allowed only for extension into 
the heads but not retraction of the filament ends back 



into the neck, we tested the stability of racquets (up to 
N = 30) to small perturbations, subject to fixed overall 
length A. We found three types of results. In the simplest 
case, the racquets are stable to any small change. These 
racquets (with N = 1,2,3,11,15,17,20...) remain ex- 
act (cr = 0). A second class is identified by stability to 
retraction but not to extension. These racquets (with 
N = 4 - 10, 13, 18, 22 ... ) develop finite (positive) over- 
hang. The remaining cases are the magic numbers start- 
ing with 12 (namely N = 12, 14, 16, 19, 21 . . .) which are 
unstable (or marginally stable) to retraction along the 
neck. A subset of these states are unstable to extension, 
and solutions with finite, positive overhang exist. How- 
ever, those cases which are stable to extension have no 
metastable solution at all. We thus conclude that no rac- 
quet solutions exist for N — 14, 16, 21, 24 . . . and these 
states are omitted from our energy spectra (Fig. ^|) and 
the series of head sizes (Fig. O). 



D. Large TV Limit and Scaling 

We perform the analogous calcu lation to that done 
for the torus states in Section VC, under the assump- 



tion that bundles form hexagonal cross sections as their 
bundling numbers N become large, to find the behavior 
of the racquet energies in the same limit. The result is 
shown as the dashed lines in Figs. |^ and [l3| To compute 
it, we assume that the large N racquet be even and with- 
out overhang (er = 0) as differences between bundles of 
nearly the same number of filaments vanish in this limit. 
This even racquet has a neck length r and a limiting head 
size Xoc for large bundling numbers n in the heads and 
N in the neck (see Fig 14). 



n = N/2 



N 



O Large N 
(even) 



FIG. 14: Schematic racquet in the limit of large N where we 
assume the symmetry of the even racquet and neglect any 
extra overhang. 

Since the size of the heads depends only on the bal- 
ance of forces at the point where the head and neck bun- 
dles meet, we can calculate the optimal head size Xoo as 
in Eq. |l3| fo r the right head of an odd racquet. As in 
Section |V Q , we determine the optimal bundling num- 
ber Af opt (A) by minimizing the energy with respect to TV, 
which yields the scaling results with prefactors as func- 
tions of A only. 



N opt « 0.303 A 4 / 5 
Xoo « 2.653 A 1/5 
Uoo « 10.482 A 3/5 



(14a) 
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FIG. 15: Spectrum of rod, racquet and torus states shown as conformational energy versus filament length in reduced units. 
Only the rod (at small A) and the tori (for all A beyond a transition point) are globally stable states. Notice the rather large 
gap between the spectrum of racquet states and the stable torus solution. Metastable torus solutions are omitted for clarity. 



Knowing the head size x<x>, we can calculate the lower 
limit of validity Ai ow in a scaling sense. This allows us to 
compare the expressions for the filament length from min- 
imization (A opt « 4.442 TV 5 / 4 ) with the length found by 
simply removing the neck altogether (Ai ow ~ 3.5757V 5 / 4 ). 
Since the optimal length A op t exceeds the minimal length 
Aiow, a large N racquet will be one with a finite neck. 
This is an important result since it hints at the evolu- 
tion of very long chains as they condense into racquets 
with increasingly larger N. In fact, we can estimate the 
growth of the neck length Too from the difference between 
the prefactors in Ai ow and A op t- Its scaling is given by 
Too ~ 0.644 A 1 / 5 . Thus, the neck grows with the same 
power of A as the heads but with a smaller prefactor. We 



may have anticipated that the growing heads provide a 
simple pathway towards the torus, as the inevitable limit 
of the heads growing at the expense of the neck. For a 
fixed filament length, the neck would have had to shrink 
to zero with increasing TV, opening the structure up to 
form a torus. For the particular racquet solutions shown 
in Fig. |we notice that the end points are relatively dense 
and represent, at times, the lowest point in the spectrum 
of states. Especially for such states, it is still true that 
their neck can shrink to very small or vanishing lengths, 
depending on A. Thermal fluctuations can then lead to 
the opening up of the neck to form a torus. Yet, even 
if the limit of large N does not provide an absolutely 
compelling pathway for the collapse to the torus, we now 
appreciate the energetics involved. 
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VII. DISCUSSION: RACQUETS VERSUS TORI 

Fig. |l5| shows the individual racquet solutions of Fig. || 
now compared to the stable torus ground states found 
in Section VA. In anticipation of these results, we de- 



scribed the lowest metastable torus state over any range 
of A as the ground state of the system. Fig. [l5| confirms 
this claim by direct comparison of racquets and tori. In 
addition, we found that the large N solutions for tori 
and racquets both grow as A 3 / 5 but with different pref- 
actors. In combination with the close agreement between 
particular solutions and the large N limit, this strongly 
suggests that the torus remains the ground state for all 
A beyond the transition point (A = 11.543). Only for 
shorter chains is the rod the ground state. 
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FIG. 16: Close-up of the rod, racquet and torus solutions in 
the region where they are closest to each other. Racquets are 
indeed never stable, though their energy is very close to both 
the rod and the tori in this region. 

There appears to be only one region where the energies 
of racquets and tori are even close, at the very low values 
of A near the transition point. Fig. |l6| shows the relevant 
region in detail. The N — 1 racquet solution comes ex- 
tremely close to the solutions for both the rod (N = 0) as 
well as the 1+ torus, but remains above. Thus the only 
stable (ground state) solutions for this system (in the ab- 
sence of thermal fluctuations) are the rod at small A and 
the tori everywhere beyond the transition point. At en- 
ergies above this ground state, we see a dense spectrum 
of metastable solutions, made up of other (metastable) 
torus (see Fig. ^ for details) and increasing numbers of 
racquet states. 

For fixed conditions we need only consider a ver- 



tical slice through the spectrum of energies. Along 
such a line, we can imagine a filament cascading down 
from an extended, rod-like configuration, through vari- 
ous metastable intermediates, while lowering its energy 
along the way. Our calculations do not of course capture 
the entire physical picture, as we neglect filament size 
in the bundling and our states are calculated in the ab- 
sence of thermal undulations. So far we have no estimate 
of the energy barriers between the metastable interme- 
diates. However, the dynamical simulation results Q 
suggest that these barriers as well as the energy gaps 
between states are large compared to fc^T: transitions 
that increase A^ are infrequent and sharp, while transi- 
tions in the opposite direction are essentially never ob- 
served. This is especially true for the transition from 
the racquet spectrum to a torus, indicating that this en- 
ergy gap is even larger for the parameters chosen in the 
simulation. This picture is consistent with the analytic 
results in Fig. [l5| which clearly shows the large gap stabi- 
lizing the ground state. The results of our analysis thus 
nicely corroborate, at least qualitatively, the results of 
our prior computer simulations as well as their relevance 
to the condensation of stiff chains. 

We would like to note that the shape of condensed fila- 
ments may depend on the nature and molecular structure 
of the condensing agent. Our study only addresses an in- 
teraction that is uniform along the filament, such as the 
effect due to a poor solvent. Other systems, with more 
point-like organizing centers, have been shown to exhibit 
intricate multi-leaf or flower patterns |32| [l7j . 

Our observations suggest that the pathway for the col- 
lapse of extended chains into condensed structures via 
intermediate racquet states is a viable, even generic al- 
ternative to the perhaps more immediately guessed direct 
winding up upon the meeting of filament ends at an ob- 
tuse angle. Some of the simulations show this latter col- 
lapse pathway, but it is much less frequent. Furthermore, 
this cascade picture through which our calculations rein- 
force and at least partially explain the simulation results, 
seems robust. We find this cascade through intermediate 
states even for a much more naive treatment of the poor 
solvent interaction used in a first pass. The individual 
curves (e.g. in Fig. |l5| ) are shifted but show a qual- 
itatively similar picture. The generic cascade through 
metastable intermediates is so dominant to be retained 
regardless of the detailed realization of the interactions. 
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